Current and universal scaling in anomalous transport 
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Anomalous transport in tilted periodic potentials is investigated within the framework of the 
fractional Fokker-Planck dynamics and the underlying continuous time random walk. The analytical 
solution for the stationary, anomalous current is obtained in closed form. We derive a universal 
scaling law for anomalous diffusion occurring in tilted periodic potentials. This scaling relation is 
corroborated with precise numerical studies covering wide parameter regimes and different shapes 
for the periodic potential, being either symmetric or ratchet-like ones. 
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Over recent years we witness an increasing interest 
in dynamical processes that occur in systems exhibiting 
anomalous diffusive behavior, possessing prominent in- 
terdisciplinary applications that range from physics and 
chemistry to biology and medicine lj . The benchmark of 
anomalous diffusion is the occurrence of a mean square 
displacement of the form (5r 2 (t)) ~ t a , where a ^ 1. 
Depending on the anomalous diffusion exponent a the 
motion can either be subdiffusive (0 < a < 1) or su- 
perdiffusive (a > 1). 

In the following we focus on the subdiffusive regime. 
Examples for subdiffusive transport are very diverse, en- 
compassing phenomena such as charge carrier transport 
in amorphous semiconductors, nuclear magnetic reso- 
nance, diffusion in percolative and porous systems, trans- 
port on fractal geometries and dynamics of a bead in 
a polymeric network, as well as protein conformational 
dynamics 0, 0, Q Another topic that recently gained 
attention is the transport of Brownian particles in the 
presence of a periodic force, relevant in Josephson junc- 
tions, rotating dipoles in external fields, superionic con- 
ductors, charge density waves, synchronization phenom- 
ena, diffusion on crystal surfaces, particle separation by 
electrophoresis, and biophysical processes such as intra- 
cellular transport 0- 0- Q 

With this work, we present intriguing results for 
anomalous diffusion and transport under the combined 
action of a periodically varying spatial force and an ex- 
ternal constant bias F. In particular, we derive a closed 
form expression in terms of two quadratures for the 
fluctuation-assisted current and its corresponding non- 
linear mobility. Furthermore, we establish a universal 
scaling relation for diffusive transport which is valid in 
tilted, corrugated nonlinear periodic potentials. 

We start out by presenting a novel derivation of 
the fractional Fokker-Planck equation (FFPE) from a 
space-continuous limit of a continuous-time random 
walk (CTRW). Our derivation involves nearest neighbors 
jumps only; moreover, it provides new insight and com- 
plements prior treatments in Refs. 

Derivation of the FFPE from the CTRW. We study 
a CTRW characterized by the probability distributions 



ipiir) for the residence times r, considering only jumps 
between nearest-neighbor sites on a one dimensional lat- 
tice {xi}, with lattice period Ax. Such a CTRW is de- 
scribed by a generalized master equation (GME) for the 
site populations Pj(t), reading 



■ equation 



{K+_ x (t - + K- +1 (t - t')P l+l {t') 



[Kt{t-t') + K-{t-t')]P t {t')}At' . 



(1) 



where the Laplace-transform of the kernel Kf(t) is re- 
lated to the Laplace-transform of the residence time dis- 
tribution (RTD) via Kf{s) = qf s^(s)/[l - ^(s)]. The 
quantities qf are the splitting probabilities to jump from 
site i to site i ± 1, obeying qf + qT = \. 

Choosing for the RTD the Mittag-Leffler distribution, 



^ i (r) = --£ a [-(. i r) a ] 1 



(2) 



one obtains Kf(s) — qfvfs 1 a . Here E a (z) = 
12^Lo 2"/r(na + 1) denotes the Mittag-Leffler function 
and the quantity v~ x is the time-scaling parameter at site 
i. The corresponding GME can be recast as a fractional 
master equation (FME) 0,^3, reading 

P l {t) = o^'i/.-i^-iW + ft+i^+iW 

- (fi+9i)Pi(t)}, (3) 

» I— a 

where the symbol oD t stands for the integro- 
differential operator of the Riemann-Liouville fractional 
derivative acting on a generic function of time as 

Him 



x(t) = 
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T(a)dt 



dt' 
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The quantities = qfvf and gi = q[ vf will be referred 
to as fractional forward and backward rates. Using the 
normalization condition for the splitting probabilities one 
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obtains that v. L = (f { +ft) 1/Q , and qf = /»/(/* + ft) and 
#r = gi/(fi + ft), in terms of the fractional rates. For an 
arbitrary potential t/(x) one can set 

/, - (K a /Ai 2 ) exp[-/5(C/ l+1/2 -Ui)], (4a) 
ffi = (K a /Az 2 ) exp[-/3(C/ l _ 1/2 - £/*)] . (4b) 

Here Ui = U(iAx) and U i±1 / 2 = U(iAx ± Ax/2), with 
U(x) the total potential, f3 — l/ksT is the inverse of 
temperature, and K(y IS the anomalous diffusion coeffi- 
cient with dimension cm 2 s~ Q . The form J3} of the frac- 
tional rates ensures that the Boltzmann relation is satis- 
fied, fc-i/gi = exp[/3(C/ l _i - Ui)]. 

By use of the Laplace-transform method one can show 
that the FME 10 can be brought into the form 

D?Pi (t)=fi-i Pi-iif) +g i+1 P l+1 (t) - (fi + ft) (t), (5) 

where the symbol £>" on the l.h.s. denotes the Caputo 
fractional derivative fill. 
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Let us introduce the finite difference operator A/ Ax, 
AP(x,t)/Ax = [P(x + Ax/2,t) - P(x - Ax/2,t)]/Ax, 
which in the limit Ax — > yields the partial derivative 
operator d/dx. Using the fractional rates (QJ the FME 
|[SJ| can now be rewritten as 



D%P{xi,t) = k c 



A:, 



(e-^-)A e ^)p (a; .^ )(6) 



where P(xi,t) = P(iAx,t) = Pi(t)/Ax. By taking the 
continuous limit in Eq. © one obtains the FFPE, 

D-P(x,t)=K a ^(e-^-^-e^P(x,t)^ , (7) 

which can be rewritten in the well-known form with the 
Riemann-Liouville fractional derivative on the r.h.s 0,0 > 



dP(x,t) _ -l-a d 



^( e -^£^P(M)).(8) 



Biased diffusion. For a constant force F the potential 
reads U(x) — —Fx and the fractional rates (0J become 
site-independent, fi = f and ft = g, satisfying the Boltz- 
mann relation, i.e., f / g = exp((3FAx) for any finite value 
of Ax. Using the Laplace-transform one finds the solu- 
tions of Eq. JSJ for the mean particle position and the 
mean square displacement of anomalous biased Brown- 
ian motion 



(x(t)) = (x(0)} + Ax(/ - g) t a /T(a + 1) 



(9a) 



(Sx 2 {t)) = (6x 2 {0)) + Ax 2 (f + g)t a /T(a + l) (9b) 
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With respect to the case of normal diffusion the expres- 
sion for the mean square displacement contains besides a 
thermal contribution oc t a also a ballistic- like term oc t 2a . 
As a consequence, a value a < 1 does not necessarily im- 
ply subdiffusive behavior. In fact, in the presence of bias 
for 1/2 < a < 1 superdiffusion takes place. 

For a finite bias F, the ballistic term in Eq. (|9b|l equals 
zero only in the case a = 1, for which normal Brownian 
motion is recovered. From Eqs. (|5J one obtains then a 
generalized nonlinear Einstein relation, which is nonlin- 
ear in force and valid for a finite space step Ax, 



{5x 2 (t)) - (Sx 2 (Q)} 
(x(t)) - (x(0)> 



= Axcoth(F/3Ax/2) . 



(10) 



In the limit F — > 0, Eq. IjlOjl yields the well-known 
Einstein relation (a = 1), n a /lJ-a(0) = between the 
thermal diffusion coefficient 



K a = T(a + 1) lim 

t — >OG 



(Sx 2 (t)) F=0 
2t a 



(11) 



and the linear mobility fi a (F = 0), with the nonlinear 
mobility fJ, a (F) = v a (F)/F being related to the anoma- 
lous current v a (see below), i.e. 



» a (F)=T(a + l) lim ^ 



(12) 



The same Einstein relation is valid also between the sub- 
mobility and the subdiffusion coefficient for any a < 1 
|l7j . as the ballistic term in the mean square displace- 
ment (|9b(l vanishes for F = 0. 

However, a relation analogous to the generalized non- 
linear Einstein relation l|10(l ceases to be valid for a < 1 
for any finite F, as the mean square displacement be- 
comes dominated by the ballistic contribution in the long 
time limit. Instead, from Eqs. 10 one obtains the follow- 
ing asymptotic scaling relation, 



lim 



(5x 2 (t)) 2r 2 (a + l) 



(x(t)Y< 



T(2a+ 1) 



(13) 



This result no longer contains the fractional transition 
rates and holds true independent of the strength of the 
bias F and the temperature T. The relation ljl3|l has been 
obtained in [18j for a continuous-time random walker that 
is exposed to a constant force. As a main finding of this 
work we prove below that this very relation holds true 
universally for the nontrivial case of tilted nonlinearly 
corrugated potentials. 

Tilted periodic potentials: fractional Fokker-Planck 
current. We next study the case of a periodic poten- 
tial with period L in the presence of a constant force F. 
Towards this goal, it is convenient to consider the FFPE 
(|T|). in which the Caputo derivative appears only on the 
l.h.s. Analogously to the case of normal Brownian motion 
0, IH ^| one obtains that the probability flux 

J a (x,t) = - KQ e-™^e^P(x,t) (14) 
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reaches asymptotically the stationary current value, i.e., 
D* a {x(t))=LJ a = v a (F) = f i a {F)F. (15) 
The anomalous current v a (F) is given in closed form by 
n a L [1 - exp(-/3FL)] 



Oa(F) =- 



(16) 



This result constitutes a first main result: It is the 
anomalous counterpart of the current known for normal 
diffusion. It obeys a form that mimics the celebrated 
Stratonovich formula for normal diffusion with a = 1 
0, fili, l20l |. From Eq. I|15|) the mean particle position 
follows as 



(x(t)) = (x{0)}+v a (F)t a /T(a+l). 



(17) 



Universal scaling in washboard potentials. We next show 
that the relation in (|13|1 is valid as well for anomalous 
transport in washboard-like potentials. We prove this by 
mapping the dynamics onto an equivalent CTRW; i.e., 
we consider a discrete state reduction of the continuous 
diffusion process x(t): To this aim, we introduce a lat- 
tice with sites {£j = jL}, located at the minima of the 
periodic part of the potential, and study the RTD tpj (r) 
for the hopping process between sites {ij}. For such a 
system, the ratio in /qj = exp(/3FL) equals that of the 
constant force case, due to the choice Ax — L. Further- 
more, the analogy between the solutions (|9a(l and Q17JI. 
both exhibiting an asymptotic power law oc t a , implies 
the same form ipj(. T ) {/T 1+a for t — > oo. In fact, for 
~ av~ a /T(l — a)r 1+a , with some suitable scal- 
ing coefficients t)j , the corresponding kernels of the GME 
obey Kp(s) = q^v° L s 1 ~ a in the limit s — > 0. Therefore, 
by making use of Tauberian theorems for the Laplace- 
transform |ll|. it follows that the asymptotic solution 
(t — > oo) is of the form JjJJl, being determined only by the 
asymptotic power law behavior of the RTD [2l| , despite 
the fact that the values of / and g depend on the chosen 
shape for the periodic potential. Because the result in 
(|13|l is independent of / and g, the scaling relation thus 
still holds true. It is universal in the sense that it holds 
independently of the detailed shape of the washboard po- 
tential, the temperature T and the bias strength F. 

Numerical verification of universal scaling. We have 
numerically tested the scaling relation (|13[) and the gener- 
alized Stratonovich formula (|16|) through the simulation 
of the fractional CTRW in tilted washboard potentials 
of various shapes and for different parameter values for 
a, F, and dimensionless T . In doing so, we not only in- 
vestigate the archetype case of a symmetric simple cosine 
potential U\(x) — cos(x), but as well a symmetric double 
hump potential U2{x) — cos(x) + cos(2x) and an asym- 
metric, ratchet-like potential U${x) = sin(a;) + sin(2z)/3. 

In Fig. we depict the scaled nonlinear mobility 
^bF ^a(F) I n a defined by Eq. (|12|l for the cosine poten- 
tial. The force is in units of the critical tilt F cr , which cor- 
responds to the disappearance of potential minima and 
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FIG. 1: (color online). The scaled nonlinear mobility 
ksT fi a {F) I ' n a is depicted for the case of a tilted cosine po- 
tential versus F/F CT . The numerics for different temperatures 
T and different a- values, varying between 0.1 — 1, fits the an- 
alytic prediction 11611 (continuous lines) within the statistical 
errors. 
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FIG. 2: (color online). Universal scaling: Asymptotic values 
of the ratio {8x 2 (t)) / {x(t)} 2 as a function of the parameter a 
for anomalous diffusion. All the points corresponding to the 
same a but different values of F and T match the function 
given in Eq. (solid line) within the statistical errors. We 
use three different temperatures: T — 0.01 with the bias F 
ranging between 0.9 — 2.0; correspondingly, T = 0.1, F = 
0.7 - 2.0 and T = 0.5, F = 0.4 - 2.0. The open triangles 
correspond to the cos-potential Ui(x), the filled squares to 
the double- hump potential U2 (x) and the open circles to the 
asymmetric ratchet potential ^3(0;), see in text. 



maxima. For a given temperature T, all values of ji a 
with a taken from the interval, 0.1 — 1, coincide with 
Eq. (|16(l (continuous lines). For tilting forces that ex- 
ceed F/F cr — 2 the dynamics approaches the behavior 
of a free CTRW being exposed to a constant bias. We 
further note that the regime of linear response at low 
temperatures is numerically not accessible. This is so 
because in this parameter regime the corresponding es- 
cape times governing the anomalous fluctuation-assisted 
transport become far too large [2p|. 

The universal scaling in tilted corrugated periodic po- 
tentials is illustrated with Fig. [21 in which the asymptotic 
ratio (5x 2 (t)) / (x(t)) 2 is plotted versus a for the three dif- 
ferent periodic potentials mentioned above. For a given 
a various data are presented, corresponding to different 
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potential shapes and values of F and T. As one can de- 
duce, these points overlap, demonstrating that the ratio 
is independent of bias and temperature, as well as the 
specific shape of U(x). At the same time, the data fit 
very well with the analytical expression l|13fl (continuous 
line). 

As detailed above, the long time behavior of the system 
is determined only by the tail of the RTD. Since we are 
interested in the asymptotic behavior (t — > oo), we have 
used in the numerical simulations the Pareto distribution 
(0 < a < 1), 

abvi _ d 1 
yjl[T) ~ (1 + bv iT y+<* ~ dr (1 + b Vi T)» ' 1 ' 

This distribution, with b = T(l — a) 1 ' 01 , has precisely 
the same asymptotic form as the Mittag-Leffier distribu- 
tion J2J • For a = 1 corresponding to a normal Brownian 
process we have employed the exponential RTD ipi{r) — 
fi exp ( — i/j t ) . In our simulations we have assumed a 
space step Ax <C L such that U"(x)Ax -C 2U'(x), ensur- 
ing the smoothness of the periodic potential. Each trajec- 
tory was assigned the same initial condition x(to) = xq. 
Then a residence time r was extracted randomly from 
the RTD (I18f) . time was increased to t\ = to +t and the 
particle was moved either to the right or left site with 
respective probabilities qf and q7 . Using this procedure 
we computed the full random trajectory of the Brownian 
particle (10 3 time steps at least). The mean displace- 
ment and the mean square displacement were obtained 



as averages over 10 trajectories. 

Conclusions. We have investigated the CTRW with 
power-law distributed residence times in a periodic po- 
tential in the presence of an external bias. The fractional 
Fokker-Planck dynamics has been derived from the cor- 
responding space- inhomogeneous CTRW in a novel man- 
ner. The celebrated Stratonovich solution for the station- 
ary current in a tilted periodic potential has been gener- 
alized to the case of anomalous transport. Moreover, we 
have proven that there exists a universal scaling law l|13|) 
— relating the mean square displacement and the mean 
particle position in washboard potentials — that does 
not involve the exact form of the periodic potential, the 
applied bias F, and the temperature T. This universal 
scaling has been verified by numerical simulations. 

Our findings for the current and the mean square fluc- 
tuations can readily be applied to the diverse physical 
situations mentioned in the introduction. The versatile 
and widespread use of the grand Stratonovich result for 
the stationary current in case of normal diffusion can thus 
readily be put to powerful use in all those multifaceted 
applications where the corresponding transport behaves 
anomalously. 
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